Topological Floquet Spectrum in Three Dimensions via a Two-Photon Resonance 
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A recent theoretical work [Nature Phys., 7, 490 (2011)] has demonstrated that external non- 
cquiHbrium perturbations may be used to convert a two-dimensional semiconductor, initially in 
a topologically trivial state, into a Floquet topological insulator. Here, we develop a non-trivial 
extension of these ideas to three-dimensional systems. In this case, we show that a two-photon 
resonance may provide the necessary twist needed to transform an initially unremarkable band 
structure into a topological Floquet spectrum. We provide both an intuitive, geometrical, picture 
of this phenomenon and also support it by an exact solution of a realistic lattice model that upon 
irradiation features single topological Dirac modes at the two-dimensional boundary of the system. 
It is shown that the surface spectrum can be controlled by choosing the polarization and frequency of 
the driving electromagnetic field. Specific experimental realizations of a three-dimensional Floquet 
topological insulator are proposed. 



Three dimensional topological insulators exhibit a va- 
riety of novel electronic properties. The most prominent 
of these are surface states, whose dispersion is that of 
massless, chiral two dimensional Dirac fermions. A dis- 
persion with an odd number of Dirac cones is unique to 
surfaces of three dimensional systems, as it is otherwise 
precluded by the Fermion doubling theorem. Such sur- 
face states were observed recently by angular resolved 
emission spectroscopy in a variety of new materials, such 
as Bia,Sbi_a; alloys, Bi2Te3,and Bi2Se3 [1-3]. The un- 
usual structure of the surface states is predicted to lead 
to unique response properties of these materials. Among 
these is the axion magnetoelectric response [4, 5], which 
arises when the Dirac cone is gapped due to breaking of 
time reversal symmetry, for example by application of a 
perpendicular magnetic field. This response is akin to 
having a fractional v = 1/2 Hall response at the surface 
of the material. Moreover, the surface states of topologi- 
cal insulators play a crucial role in proposals for quantum 
interference devices. Most notable of these is the possi- 
bility to realize and manipulate Majorana Fermions [6], 
which have important applications for topological quan- 
tum computing. 

The topological behavior of electrons is emerging as a 
promising resource, and therefore it is imperative that we 
understand all ways to induce it. In this manuscript we 
explore the possibility of dynamically inducing a three di- 
mensional topological spectrum, surface states included, 
starting with a trivial (non-topological) bulk insulator. 
The idea of inducing topological order with periodic mod- 
ulations of a Hamiltonian was explored in Refs. 7-10. 
The Floquet spectrum of a periodically driven system 
was shown to exhibit a variety of topological phases, in- 
cluding one that exhibits a single Dirac cone [7]. Phys- 
ical examples include graphene which is expected to ex- 
hibit a quantum Hall effect when subjected to radia- 



tion [9, 11, 12], and spin-orbit coupled semiconductor 
heterostructure (such as HgTe/CdTe wells), which can 
be turned from trivial to topological using microwave- 
teraHertz radiation [10], and vice versa [13]. 

In this manuscript we demonstrate how a "time rever- 
sal invariant" three dimensional topological spectrum can 
be induced in trivial insulators using electro-magnetic ra- 
diation. This problem is a non-trivial generalization of its 
2D analog [10]. Roughly speaking, a topologically trivial 
band structure can be turned topological by mixing the 
valence and conduction bands, for instance by radiative 
transitions. In 3D, the radiation has to be carefully tai- 
lored such that it produces a non-vanishing band inver- 
sion matrix element in a closed 2D surface in momentum 
space, and, as we shall see, must obey additional topo- 
logical and symmetry constraints. The polarization and 
frequency of the driving electromagnetic field allow for a 
detailed engineering of the surface states, including the 
possibility to carefully tune a gap in the Dirac cone at 
the surface. 

First, let us develop our ideas within a simple generic 
band structure. We consider an effective low energy 
model near the F (k = 0) point [3]. The four states 
near the Fermi energy at the F point are denoted using 
two quantum numbers, corresponding to spin a =t, i and 
parity r = +, ^. Time reversal symmetry is represented 
by T = ia-ylC, where IC is the complex conjugation opera- 
tion. Inversion symmetry is represented using I = I^Tz- 
We study a Hamiltonian of the form 

H = 5(k) • 7 (1) 

where 7 = (71,72)73; 75) are four Dirac matrices, which 
we represent by 7^ = CTj (g) with i = 1, 2, 3, 75 = / (g) r^. 
The remaining dirac matrices are defined as 74 = 7 (g) Tj; 
and 7ij = — 2i[7i,7j]. In Eq. (1) and below, we denote 
3-dimensional (space coordinate) vectors, such as the mo- 
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mentum, k, in bold symbols, while the 4-diniensional vec- 
tors are denoted with a vector symbol D, and D for unit 
vectors. Writing D{k) — (d(k); £>5(k)), we note that 
the Hamiltonian (1) has both time reversal and space- 
inversion symmetries under the restriction that the vec- 
tor d(k) be odd under inversion, while I?5(k) is an even 
function. 

Time- reversal- invariant (T^ = — 1 ) band insulators 
in three dimensions admit a Z2 classification [14-16], 
falling into two categories - either topological, or trivial. 
The model in Eq. (1) can describe both phases, depend- 
ing on the choice of parameters. Let us emphasize that 
while Eq. (1) does not describe the most general Hamil- 
tonian with time reversal symmetry in three dimensions, 
this effective model does span a wide variety of realistic 
systems and allows for a relatively simple visualization of 
the Z2 topological invariant that we now focus on. 

At each momentum, k, the spectrum of the Hamil- 
tonian (1) is doubly degenerate. Its eigenstates Vk 
of (1) are also the eigenstates of the rank-two projec- 
tors P±(k) — ^[I zL -D(k) • 7] onto the valence (— ) and 
conduction bands We can parameterize the unit vec- 
tor, I?(k) = £'(k)/|£'(k)| (lying on a three-dimensional 
sphere, S^), using two polar angles, 6 and ^, and an axial 
angle, (/>. We define 6 as cot(6»u) = i:'5(k)/|d(k)|. The 
angles ^, correspond to the spin direction, by the unit 
vector d(^k, 0k) = d(k)/|d(k)|. Note that there exists no 
global coordinate map on and in our case, d(^k, 0k) 
remains undefined at = 0, tt. 

Using the above parametrization, the topological in- 
variant for Hamiltonians of the form (1) can be calculated 
by considering the map from the three dimensional Bril- 
louin zone (BZ), which is a three dimensional torus, T'^, 
to S^. The only topological invariant of this map is an 
integer which counts the number of times the map wraps 
the target space S^, also called the degree of the map. 
Two Hamiltonians of the form (1), for which the degree 
of the map differs by 2, can be deformed into each other 
without closing the gap in the spectrum [17], by adding 
terms which digress from the form (1). Therefore, the Z2 
classification of the insulator (1) is given by the degree of 
the map mod 2, with an even degree corresponding to a 
trivial insulator, and an odd degree to a topological one. 

Near the F point, an expansion to order k^ yields 

D{k) = {Ak;M - Bk^), (2) 

where spherical symmetry can be assumed for conve- 
nience. The topological phase of the Hamiltonian (1),(2) 
occurs for M/B > 0. In this case, the angle, 6, changes 
from 61 = at k = to 6* = vr at [k[ > ^JM/B. For each 

< < TT, the vector d wraps a sphere of two dimen- 
sional unit vectors, S"^ , which corresponds to a "latitude" 
on the target space S^. Therefore, the target space of 
the map is covered once in the topological phase. 

Consider now the case of M/B < 0. In this phase, the 
valence band is characterized by < 6* < 7r/2, while the 
conduction band corresponds to Tr/2 > 6 > n. Therefore, 








FIG. 1. The two paraboloids represent the dispersion relation 
e±(k) of the valence and conduction bands for the Hamilto- 
nian (1) in the trivial phase, projected on = 0. Each energy 
is doubly degenerate. At each momenta, the eigenstates of (1) 
can be represented by their direction in r space, tan{6) = 
\{r^)\/{Tz), and spin direction d(C, 0) = (-1)"'^"^-) (o-) . The 
color scheme of the paraboloids represents the value of the 
angle 6. The spin direction on spheres in momentum space, 
|k| = ko is depicted on the right. A two photon resonance 
connecting the valence and conduction bands occurs at such 
a sphere in momentum space, and is represented by the green 
circles on the parabolas. 

the degree of the map from the BZ to is zero, lead- 
ing to a trivial insulator. We note that terms involving 
other 7 matrices can be added to H{\C) in Eq. (1), while 
keeping time reversal symmetry. Then, the Z2 topologi- 
cal invariant cannot be calculated using the above simple 
considerations. However, as long as the added terms do 
not cause the gap in the spectrum of i/(k) to close, the 
Z2 invariant does not change. 

Having reviewed the construction of the Z2 topologi- 
cal index within the effective model (1), we now discuss 
a non-equilibrium case, where a time dependent pertur- 
bation is added to Eq. (1) initially in a trivial phase. As 
we shall see below, to induce a topological band struc- 
ture in three dimensions using electromagnetic radiation, 
the external field has to satisfy non-trivial constraints. 
We first illustrate these requirements and elucidate the 
physics in the framework of the effective model with a 
generic time-dependent perturbation of the form: 

V{t) = Re (v^e*"*) • 7 (3) 

where I^ is a complex, fixed, four-component vec- 
tor. According to Floquet theory, the unitary opera- 
tor describing the time-evolution of the system, U (t) ~ 

V exp ^— i Jj*^ dtH{t)j , can be written as 

U{t) = W{t) exp [~iHp(]s.)t] (4) 

where W is a unitary matrix satisfying W{t + T) = 
W{t), T = 2tt/uj and i7i?(k) is a time-independent Flo- 
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quet operator. The cigcnstates of Hp are solutions of 
[-idt + H{t)] = of the form = e-*'^*$(t), where 
$(f) is periodic with the period, T, and e are called the 
quasi-energies. In the following, we shall take the time 
independent piece in the Hamiltonian, i?(k), to be in the 
trivial phase, and show that a topological spectrum for 
the Floquet operator iJ^(k) can be achieved nonetheless, 
by choosing V{t) appropriately. 

In three dimensions, and in the absence of particle- 
hole or sublattice symmetries, a topological spectrum 
requires time reversal symmetry [14-16]. The instan- 
taneous Hamiltonian _H"(k,i) = -ff(k) + V{t) docs not 
necessarily posses such symmetry. However, if we can 
satisfy 



TH{\<i,t)T-^ = H{\<i,-t + T) 



(5) 



for some fixed but arbitrary parameter, r, then the Flo- 
quet operator HpCk) is invariant under an effective time 
reversal symmetry [7, 10]. The time-reversal constraint 
in Eq. (5) is satisfied if arg(Vi^2,3) = arg(V5) -I- tt. 

The effect of the time-dependent potential becomes 
apparent in rotating frame, where the lower band is 
shifted up by lo. Wc transform the Hamiltonian to 
such a rotating frame via the unitary transformation 
C/(k,t) = P+(k) -I- P_(k)e-*"* (here P±(k) are the pro- 
jectors defined above). The resulting interaction-picture 
Hamiltonian reads 



Hi(k,t) = e+{k)P+{k) + [e. 
+ U{t)V{t)UHt), 



(k)+a;]P_(k) 



(6) 



where e±(k) = ±|D(k)| are the band dispersion relations 
of ff(k). The Floquet operators corresponding to Hj{t) 
and to -ff(k) are identical, up to a (time-independent) 
unitary transformation. 

Consider the first two terms in Eq. (6). In the inter- 
action picture, the two bands intersect on a two-sphere 
in the BZ, |k| = feg, if the driving frequency to is larger 
than the band gap, 2M. We denote this two-sphere by 
S, and we depict it in Fig 1. The driving term (third 
term in Eq. (6)) opens a gap in the quasi energy spec- 
trum at ]k| = ko. For momenta away from S, the effect 
of the driving terms in Eq. (6) can be neglected [18], and 
the quasi-energy states arc roughly cigcnstates of P± (k) . 
Notice, however, that the bands of Hj are inverted: the 
projector P{ (k) onto the lower quasi-energy band of Hj, 
corresponds to P+(k) near the F point, and to P-(k) 
for |k| » ko (and vice versa for P_[^(k)). The projectors 
P^(k) smoothly interpolate between these two limits. In 
order to find the degree of the map from the BZ to 
that the projectors Pj-(k) define, we need to examine 
their properties on the resonance sphere S. 

For values of k near the resonance sphere, the rotat- 
ing wave approximation (RWA) reveals the effect of the 
driving field. In this approximation, only time indepen- 
dent terms are kept in the last term of Eq. (6). This 
yields 



Vrwa 



^(^P+VP- + P-VKP+y V = V-i (7) 



The topological properties of Hi are intimately re- 
lated to the transformation properties of Vrwa under 
the group of spatial rotations in three dimensions. We 
shall take the action of this group in spin space to be 
exp(im • <t), where m parameterizes the axis and angle 
of rotation. Note that while 75 transforms trivially under 
spatial rotations, (71, 72- 73) transform as a vector. Since 
to second order in k the projectors P±(k) are scalars, the 
transformation properties of Vrwa at this order arc de- 
termined by those of V. Therefore, to second order in k, 
Vrwa contains both scalar and vector representations of 
the rotation group. As wc explain below, a scalar Vrwa 
yields Pi-(k) with a non- trivial topological Z2 invariant, 
as it maps the resonance sphere <S to cover each spin di- 
rection once. However, a purely vector contribution to 
Vrwa yields a Pj. (k) which is topologically trivial, as the 
map from S covers only a partial cap of spin-directions. 

The explicit form of Vrwa is 

Vrwa = 1(^-7 + AIni{ V,- Hj) (8) 



where 



V± = Re{V} - (Re{V} • D)D. (9) 



Two illuminating cases are: (i) V = and (ii) V = 

ViXi, with Vi and V5 real. This gauge choice allows us, 
using Eq. (9), to approximate P±(k) « ^(1 ± n(k) • 7), 
and study the map from the BZ to defined by n(k). 

In case (i), V is a scalar under spatial rotations, and 
hence so is Vrwa = V_l • 7. Therefore, the first three 
components of V± must be proportional to k while the 
fifth components is a scalar. Indeed, using Eq. (9) we get 



= (- 



£)2 



P»2 



(10) 



where D — |P)(k)|. On the resonance sphere S, we have 

n{\c,t) = V_l/|V_l|, which maps <S to an S"^ sphere with 

s 

fixed latitude 6{ko) in the target space S^. Therefore, 
n(k) defines a map from the BZ to S^, which maps the 
r point to the south pole, 5*^ spheres in the BZ to 
spheres ("latitudes") on S^, and maps k ^ ko towards 
the north pole of S^. It is therefore a map of degree one, 
which implies a topological spectrum. 

In Case (ii), however, V and Vrwa give a vector repre- 
sentation of spatial rotations. Therefore, the first three 
components of V± must either be scalars or belong to 
a spin-2 representation of spatial rotations, while the 
fourth and fifth components must belong to a vector rep- 
resentation. Indeed, an explicit calculation gives 



V± = Vix- 



A^kxk Ak, 



£>2 



D D 



(11) 



Clearly, on S, the vector V± does not wrap around an S'^ 
sphere on S^, and the resulting map defined by fi(k) is 
topologically trivial. 
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More generally, Vrwa will be a superposition of a 
scalar and vector component. Whether a topological 
spectrum is obtained depends on the relative magnitude 
of the two components. 

Let us now apply the geometrical considerations lead- 
ing to a topological Floquet spectrum in 3D, to oscillat- 
ing electromagnetic fields. The electric field operator is a 
vector under spatial rotations. Therefore, in light of the 
above discussion, it first seems impossible that it could in- 
duce a topological Floquet spectrum. This can be reme- 
died by considering multipole tensors of the electric field. 
The quadrupole tensor EiEj can be decomposed into a 
scalar, which is given by its trace, and a tensor giving a 
spin-2 representation of spatial rotations. In the follow- 
ing, we show how to use the scalar part of the quadrupole 
tensor in order to induce a topological Floquet spectrum. 
The scheme involves choosing a frequency to which sat- 
isfies M/2 < CO < M. Such a choice for the frequency 
precludes a resonance induced by a single-photon tran- 
sition, which is linear in the driving electric field E, but 
allows for a two-photon transition, which is second order 
in E, and therefore involves the quadrupole tensor. 

Note that in order to satisfy the time reversal sym- 
metry constraint, as defined by Eq. (5), the oscillating 
field must be linearly polarized [10]. An ellipticity in the 
polarization of the driving field leads to interesting ef- 
fects, which will be discussed later. We choose a gauge 
E = dtA, (/) 0, whereby the Hamiltonian becomes 
H = D{\s. + A{t)) ■ 7. Choosing A = Aq cos{ujt)-k, we 
obtain the Hamiltonian 

H{t) = i5'(k) • 7 + • 7C0s(wt) + V2 ■ 7COs(2wt), (12) 

with Vi = AAqx- 2BAokxX5, V2 = -^BAlii^, and 
i5'(k) = [Ml; M - BAl - Sk2^ . Below we sketch the 



outline of the calculation, the details of which are given 
in Appendixes A and B. 

The two-photon resonance is a second order process 
in the electric field. In the chosen gauge, a 2w-term arises 
both directly from V2, and from a second order process 
in Vi. The contribution of both terms to Vrwa scales, 
to lowest order in the light intensity, as Aq/M. In or- 
der to calculate their effect, we perform two consecutive 
unitary transformations (see Appendix B for an alterna- 
tive derivation). We first perform a transformation to 
a frame rotating with frequency co, of the form leading 
to Eq. (6). The resulting interaction picture Hamilto- 
nian does not contain any resonances. We diagonalize 
its time independent terms, which yields new eigenvalues 
and projection operators, which we denote by (k) and 
P^\\s.) respectively. A second unitary transformation. 



?72(k, t) = P^'' + P^'^' exp(— iwt) yields a new interaction 
picture Hamiltonian 

H2 = eV'Hk)pW + (£L'^(k)+u;)pW+f/2V(^)(0t^2 (13) 

where V^^^{t) denotes the time dependent terms resulting 
from the first transformation. 



,(1) 



After the second transformations, the two bands cross 
at a surface S with £+^(k) = (eL^^(k) -|-w) and the topol- 
ogy of a sphere. We now employ the rotating wave ap- 
proximation, which yields on the resonance surface. 



-n2,RWA|5 — 2^± 



(14) 



The vector is defined, to lowest order in ^0 and 
k, as in Eq. (9) with the replacement D D^^^ and 
V = (^{Vi - V2) ■ D + V2 (see Appendix A 

for details). The contribution to v}^^ which allows for a 
topological map comes from the term corresponding to V2 
in Eq. (12). A quick inspection shows that l^-f is a scalar 
under spatial rotation, and after the two transformations 
yields a contribution to vj^^ of the form appearing in 
Eq. (10), with corrections to its spatial components of 
higher order in k and ^o- The contribution to vj^^ com- 
ing from Vi in Eq. (12) only leads to an anisotropy of the 
gap in the Floquet spectrum, and does not change the 
topological properties of V^^ on the resonance sphere <S. 

The first three components of the vector field vj^^ 
are plotted in Fig. 2. Clearly, they map <S to a single 
covering of the unit sphere. Note that the 5*'' component 
of vj^"^ is not constant on S. This does not change the 
degree of the map from the BZ to S^, and the Floquet 
spectrum is characterized by a non-trivial Z2 invariant. 
The magnitude of the gap on the resonance surface is not 
isotropic. This is expected, as choosing the polarization 
direction breaks the rotational symmetry of the problem. 
On the resonance, the Floquet spectrum is fully gapped, 
where the smallest gap occurs for = 0, and is given by 



gap 



ABA^ 
4M 



(15) 



where k^ is the resonance wave vector. 

One of the most striking consequences of the topolog- 
ical band structure in three dimensions are the mid-gap 
surface modes , which are characterized by a single Dirac 
cone [19]. Likewise, a striking consequence of the above 
considerations are the appearance of surface modes in the 
presence of the driving electric field. The surface quasi- 
energy states appear inside the Floquet quasi-energy gap 
and are characterized by a single Dirac cone. 

To demonstrate this, we use exact numerical methods 
to study the Floquet problem of a tight binding model 
corresponding to the Hamiltonian of Eq. (1). We con- 
sider a finite slab with vanishing boundary conditions at 
z = 0,L and a driving electric field polarized in the x 
direction. The quasi energy and momenta in the x and y 
directions are good quantum numbers. In Fig. 3, we plot 
the quasi energy spectrum inside the quasi-energy gap, 
which clearly exhibits a single dirac cone. Note that the 
Dirac cone is not isotropic, resulting from the choice of 
the polarization along the x direction. 

Broken time reversal symmetry on the surface of a 
three dimensional topological insulator leads to a gap at 
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FIG. 2. The first three components of the vector field l7[^^(k) 
resulting from a two photon resonance of a linearly polarized 
electro-magnetic field. The magnitude of the plotted vector 
field, on the resonance sphere S (depicted), gives the gap 
in the Floquet spectrum. Its direction determines the spin 
direction of the quasi-energy states on 5. The map from 5 
to the two-sphere 5*^, defined by these three components, is 
of degree one. The resulting gap in the Floquet spectrum is 
not isotropic, due to the necessary choice of the direction of 
polarization of the electric field. 



the Dirac node of the surface states. Such a gap en- 
tails unique consequences in terms of transport proper- 
ties of the surface states, as it leads to a half integer 
quantum hall effect on the surface. Remarkably, a gap 
at the Dirac node of the surface Floquet spectrum can 
be easily controlled by an appropriate tuning of the po- 
larization of the electromagnetic driving field. Consider 
an elliptically polarized electric field, E = ReEe"^' with 
£ — E/|E| = (x-|- iSy). The electric field operator is 
even under time reversal. Therefore Eq. (5) cannot be 
satisfied for any 6^0. The gap in the Dirac cone on the 
surface depends quadratically on the elipticity parameter 
S. 

The analysis above demonstrates how to obtain a 
"time reversal" invariant topological spectrum in three 
dimensions, by periodically driving a topologically triv- 
ial system. Our approach could be applied in a va- 
riety of quantum systems, e.g., cold atoms with syn- 
thetic spin orbit couplings [20, 21]. The most natural 
experimental application of our results are in semicon- 
ductors with appropriate properties, namely, a spectrum 
with a narrow direct band-gap occurring in one (or an 
odd number) of points in the Brillioun zone. Candidate 
materials are Sb2Se3 in the rhombohedral crystal struc- 
ture [3], GeSb2Te4 [22]. Moreover, Heusler compounds 
[23, 24] with applied strain exhibit four bands and a nar- 
row bandgap near the T point, and are therefore excel- 
lent candidates for our proposal. In the materials men- 
tioned above, a sizable gap in the Floquet spectrum, on 
the order of lOK, can be achieved using electric fields of 
lO'^V/m. 




FIG. 3. Quasi-energy spectrum of the Floquet operator corre- 
sponding to Eq. (12), in the slab geometry: periodic boundary 
conditions in the x, y directions and vanishing ones in the z 
direction. In a) and b), we plot the quasi-energy spectrum in- 
side the gap of the Floquet spectrum, as a function kx,ky, for 
different polarizations of the driving field, a) Linearly polar- 
ized electric field, e = x, yielding a Dirac cone on the surface, 
b) A small gap in the surface spectrum is opened for an el- 
liptical polarization with small ellipticity, e = x + i5y, with 
5 = 0.05. Note that in both a) and b), the Dirac cone is not 
isotropic, c), d) Quasi-energy spectrum, showing bulk bands 
and edge states, c) Spectrum as function of k^, for ky = 0, 
with e = X. d) Quasi-energy spectrum as a function of ky, 
for kx = 0, e = X -I- iSy. The inset magnifies the gapped 
Dirac surface spectrum. Note that the the gap in the Floquet 
spectrum is enhanced for kx > 0, as evident by comparing the 
spectra in c) and d) 



The surface modes in the Floquet spectrum predicted 
in this work could be directly probed using photoemis- 
sion spectroscopy [25, 26]. Moreover, a gapped surface 
Floquet spectrum can be detected using the Kerr and 
Faraday effect [27], whereas in the material these effects 
should be negligible absent the driving. The emergent 
Floquet spectrum discussed above is that of a band in- 
verted semiconductor [28] , and in any steady state we can 
expect a finite density of particles and holes in the Flo- 
quet bands [29] , [30] . One concern is that these free charge 
carriers may attenuate the driving electromagnetic field. 
For a lower bound of the attenuation length C we assume 
that all electrons participating in the band inversion act 
as free charge carriers. This leads to C ^ ^ ^ 1/ ^JJm, 
with carrier density n ^ fcg and mobility /i. Since the 
localization length A of the surface modes scales linearly 
with k^^ , a parameter regime in which A < ^ could be 
found, according to material parameters. 
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Appendix A: Topological properties of the 
two-photon resonance 

In this section we study the topological properties 
of the Floquct operator corresponding to an insulator 
driven with an electromagnetic field whose frequency al- 
lows only for a two photon resonance, M < cj < 2M, 
where 2M is the bandgap of the insulator. The two- 
photon resonance is a second order process in the electric 
field. In the chosen gauge, a 2u;-term arises both directly 
from V2, and from a second order process in Vi. Both 
terms therefore yield contributions of order Aq. In or- 
der to consider both terms on an equal footing we shall 
perform two consecutive time dependent unitary trans- 
formations, where each transformation is characterized 
by the frequency w. In order to analyze the resulting 
time dependent Hamiltonian, we shall employ the rotat- 
ing wave approximation and expand to lowest orders in 
k and ^o- 

We first introduce some useful notations. We decom- 
pose any four dimensional vector V into the components 
parallel and perpendicular to a four dimensional unit vec- 
tor D as 

V^n = V-{v-D)D, (Al) 

and 

V^l^ = {v-d) D. (A2) 

For notational convenience, we shall relabel -D'(k) — >■ 
D{k), c.f. Eq. (12). The first rotating wave transforma- 
tion is done via the unitary 

[/i(k,t) = P+(k) + P_(k)e-'"*, (A3) 

which leads to the Hamiltonian in the first rotating frame 
given by 

H, = Ui{\^,t)H{t)Ul{]^,t) 

= (5(k) + IX?^^) . 7 + c^P-(k) + V(i)(t).(A4) 

In the above, the time dependent part is given by 

V^'Kt)= (yi||^ + il72i^).7Cos(c.i) 

+ hiV2±D)Alij M^t) + Vit). (A5) 



Note that V*^"^^ (t) contains terms with frequencies cj (the 
first two term in Eq. (A5)) and 2ui, 3w (corresponding to 
V{t), the third term above). 
It is convenient to define 

^(i)(k) = Ae(k)I)(k) + ifi^^(k), (A6) 

with 

Ae(k) = 6(k) - ic, e(k) = |^(k)|, (A7) 

which enables us to write Eq. (A4) as 

Hi=D^^\k)-j + V^^\t) (A8) 

The time-independent part of Eq. (A4) can be ex- 
pressed using eigenvalues £±^(k) and projection opera- 
tors P^±\k). We now perform a second rotating wave 

transformation 

U2{k,t) = Pl^\-k) + P^^\k)exp{-iojt), (A9) 
which yields the Hamiltonian in the 2'^'* frame, 
H2 = e^'^(k)Pj.')(k) + (eW(k) + w)P«(k) 

+ c/2(k,t)((y,||^ + iK,^^) • j)ul{k,t)cos{ut) 

+ U2{k,t)(^{lV^^j^),I),j,,)uUKt)sm{ivt). (AlO) 

In the above, we have omitted from H2 the term 
U2V(t)ul which does not contain any time independent 
contributions to H2, and therefore does not contribute to 
the two photon resonance. 

After the second transformations, the two bands cross 
at a surface S for which (k) = (eL^^ (k) -I- to) . We now 
employ the rotating wave approximation. The contribu- 
tion coming from the second term in Eq. (AlO) can be 
deduced from inspecting Eq. (8) and (9) in the main text. 
The contribution arising from the third term in Eq. (AlO) 

yields a term of the form ^Di{V2j_jj)jDk[jij,'^k]- Some 
algebra reveals that to lowest order in ^0 and k, the 
two terms in Eq. (AlO) involving T^2_Lri yield the same 
contribution to the rotating wave approximation. 
Therefore, on the surface <S we have, 

H2,KWA\s = lv[%,r^ (All) 

where the vector t^^^fij is defined as in Eq. (Al) by re- 
placing 

D^D^'\ ^^l^^''=V^i||D + W (A12) 

In order to achieve a topological Floquet spectrum, 
the vector field ^^(ij needs to map the resonance surface 

S to an S"^ sphere on the three dimensional sphere S^. In 
the following we shall show that this is indeed the case. 
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As a first step, we inspect the contributions to V'^^\ 
which arise after the first unitary transformation. Keep- 
ing only terms only up to second order in k, we find 



Au{A- - 2DM) 



IIID 



M 



and 



V --^(k 



M 



(A13) 



The vector field V21. clearly maps a sphere in the BZ to 
an 52 sphere on the target space S^. This is of course 
expected from noting that V2 ■ 7 is a scalar under spatial 
rotations. However, we are interested in its contribution 
to the two-photon resonance, i.e., to Eq. (All). 

To this end, we note that V^^f) is orthogonal to D by 
construction, and D^^) =b + o(^) w £» + o(J^). 



Therefore, the correction to 



into the expression for ^^^(i), are of higher order in ^o- 
Explicitly, we have 



2±D' 



when it is inserted 



V - -IJ^i) = 



— — V21. ■ Vi 
2|£)(i)| 

1 ABAl ( 



1 A^BAl 



A 



AB ,\ 
2— k^ fc^ 
M J 



Ae(k) 4M 

where we have kept terms only up to order k^ and Aq. 
Therefore, the final contribution of V2 to ^^^(i) to this 
order is 



2±D 



( 



A^BAl 



V4AfAe(k) 



(A14) 



The second term in the above equation correspond 
to a correction to the spatial (1 — 3) parts of ^2_Lr)' 
Eq. (A13), of order \s? and ^q. The spatial part of l^iij 
are originally of order k and A^. Therefore, to lowest 
order in k and Aq, this correction does not alter the 
topological property of which we shall describe 

below. 

We now calculate the contribution of Vi to ^^^(i, 7 
which will turn out to be of the same order as the con- 
tribution of V2, c.f. Eq. (A14). 

First, we note that 

V - ■ - — ' 



\DW\ 



(A15) 



From Eq. (A15), and for Ae(k) > |V^j^£)|, we have 



2Ae 



Using the definition of Vin, Eq. (A2), we see that the 



total contribution of Vi to "^^^^(1) , to order Aq is 



yiii-(yii|.^«)i>w 



2Ae 



i± 



(A17) 



Therefore, this contribution to 'l^^^(i) is of the same or- 
der in the driving field as the contribution coming from 
^2_L£)' (^13) and its inclusion was necessary for com- 
pleteness. Note that Vi± = AAox + o{\\i\), and therefore 



2Ae 



1±D 



Al{A^-2BM) 
2MAe 



-Ak, + o{Al,\k\^) 
(A18) 

while the y and z components of Eq. (A17) are of order 
Aq and k'"*. Note that the xr, component of Eq. (A17) is 
of order Aq and k^. From Eq. (AI6), we sec that the Vi 
term also contributes terms of order Aq and k^ to ^^(i) • 
All of the above higher order corrections do not change 
the topological properties of to lowest order in 

Ao and k. 

Summing up both contributions to l^^^(i) ) we have, 
to lowest order in k and 



= ( 



Vi±) 



2Ae 

2MV--(^^- 
while the other two spatial components are given by 



^1(b. 



2BM)/AejAk:,iA19) 



^Ak 
2M 



a = y,z 



(A20) 



From Eqs. (A19,A20), we sec that the vector field 
^±£)W maps the resonance surface S to an 2-sphere on 
the target space S^. This 2-sphere is not at a constant 
"latitude", i.e. its 6 coordinate on is not constant. 
Importantly however, this 2-sphere winds around the 
poles of S^, i.e., it is an incontractible sphere on the 
space S"^ \{N Li S), the space of with the north and 
south pole removed. 

From the above considerations, we see that i?2(i), 
given in Eq. (AlO) can be characterized by projection 

operators of the form -^(k) = |(1 ± ^2(k) • 7). The 
unit vector n2(k) defines a map from the BZ to with 
the properties: (i) For regions in the BZ near the F point, 
fi2(k) Ri —D{k), and therefore it maps 2-sphercs in the 
BZ to 2-spheres on 5^, close to its south pole; (ii) Maps 
the sphere S in the BZ to an S"^ sphere on which 
is incontractible on S*"^ \ (A^ U S) (as discussed above) ; 
(iii) For large values of k, n(k) w -D(k), therefore these 
arc mapped to 2-spheres close to the north pole of S^. 
From continuity of n2(k), it must therefore define a map 
of degree one. This implies that the Floquet operator 
corresponding to H2{t) and H{t), has a non trivial Z2 
topological invariant. 
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Appendix B: Virtual absorption perturbation theory 



The two consecutive RW transformations are very reminiscent of a perturbation expansion. One RW transforma- 
tion fails to produce a dcgcnracy, and therefore a second transformation is necessary to expose the role of two-photon 
processes. In this section we will show how indeed such processes can be analyzed as virtual absorption processes, 
and derive a formula which replaces second-order degenerate perturbation expansions. 

The first step involves mapping the time-dependent Floquet problem to a time-independent problem using an 
auxiliery degree of freedom. Let us introduce an additional Hilbert space which serves as a counter of photons 
absorbed (for the experts, we note that this auxiliery variable is just a way of keeping track of the Floquet block 
index). We introduce an infinite lattice for a single particle, which we denote F, with states \n)p; n is essentially 
counting the number of photons absorbed by the system. The original system has a Hamiltonian which is split to 
time independent 'Hgys and time dependent pieces: 



H{t) = Hsys + 6e''^' + 6^e- 



(Bl) 



we now replace the time dependent terms with hopping terms for the register particle F. We also add a diagonal 
energy term that determines the energy of the F states. The Hilbert space after this mapping is a tensor product 
state between the F-states and the system's states. The time dependent Hamiltonian is therefore replaced by an 
operator that acts on the larger Hilbert space, 



y-F = y-sys + ^ (O |n + 1)^ p {n\ + 6^ \n)p p{n+ 1|) + H^:, 



(B2) 



with 



Hu = ^noj \n)p p (n| . 



(B3) 



and "Hsys = F^ys SiIf- 

To retrieve the original Hamiltonian, Eq. (Bl) all that is necessary is to initiate the auxiliary F-states in the 

zero-momentum state: 



(B4) 



with N providing a normalization. 

Our first claim is that the time- independent Hamiltonian, Eq. (B2), initiated with \'ip)p has the same propagator 
for the system as the one for the original Hamiltonian, Eq. (Bl): 



U{t) = V 



exp 



i / dtH(t) 
Jo 



-F (V"! exp[i'H(^f] exp {-iHpt] 



IF ■ 



(B5) 



where V denotes path ordering. To show this, we first move to the interaction picture in terms of the states F. More 
precisely we consider: 



We note that 



U{t) = e*«"* • e-'«-*. 



Since T-Lf — "H^; = T-Lsys + Hqf, we write: 



€^'^-\Hf - n^)e-'*^-' = Hsys + HoF{t) 



(B6) 



(B7) 



(B8) 



with: 



noFit) = e'^-'Yl \n + l)p p (n| + 6^ \n)p p{n+ 1|) e"''^-* = ^ (o |n + 1)^ F {n\ e^"* + 6^ \n)^ p{n+l\ e"' 

n n 

(B9) 
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U{t) is easily solved to be: 



exp 


h/ 




. Jo 



dt{U,y,+U0F{t)) 



And therefore, the identity (B5), which we are trying to prove, becomes 

U{t) =F {i>\U{t)\i,)^. 



(BIO) 



(Bll) 



Now that we have essentially eliminated Tiu, from the expression for U{t), the only operators relating to the F-states 
remaining are the hopping operators |n + 1)^, p {n\ and Yl \''^)f f ('^ + ^l- These operators are simple to handle 

n n 

since \'4')p is an eigenstate of both, with eigenvalue 1: 



^ |n + l)p p {n\ iI})f = ^\n)p p{n+l\ 'iP)f = \'4')f- 



(B12) 



Thus we can write: 



and also: 



noFit) I V)f = Wf (Oe'"* + Ote--*) = \^|;)^ Ho{t) 



U{t)\ 



exp 



I dt{Hsys+H0F{t)) 

Jo 

i I dt{Hsys+Ho{t)) 

Jo 



\i>)F 



'U{t) 



(B13) 



(B14) 



which confirms Eq. (Bll), and therefore completes the proof of the mapping. 

To conclude this section, we note on the correspondence between the formalism presented above and the Floquet 
theorem 



U{t) = W{t) exp [-iHFt] 



(B15) 



where W{t + T) = W{t) and Hf is an operator acting on the system Hilbert space only, see main text, Eq. (4). The 
correspondence is given by noting that choosing W{t = 0) = I gives 



exp [-iHpt] = (V'lir exp [-iUFt] \il))p 



(B16) 



1. Elimination of single photon processes 



The auxilicry F-statcs formalism allows accounting for a sequence of virtual photon absorptions, by systematically 
eliminating the F-states associated with intermediate parts of the process. In the case we considered, for instance, 
there is no resonant single photon process. Therefore, if we start the system and F-state wave function with only an 
even number of photons, odd-photon F states will only appear with a STipprcssed amplitude since they have a large 
energy mismatch with the initial states of the wave function - they must be about an w away. 

Accounting for two-photon processes in our system is thus possible along the lines of ordinary second-order 
perturbation theory. We start by consdiering the propagator applied to the low-energy subspace, and read-off the 
effective hamiltonian that emerges after resumming connected diagrams. In our case, the low-energy subspace of the 
F-states is the superposition of all even states: 

|V'e.en)^ = ^^|2n)^. (B17) 
n 

The effective two-photon propagator is then given by: 

U{t) » C/2(t) =F {^even\U{t) \^even) F ■ (B18) 

Next, we need to expand the interaction Hamiltonian in powers of T-LoFit). 
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Before carrying out the expansion, let us move to the interaction picture of Hsys as well: 

ei-Hsyst^^f^t) «^ (Ve^enl ^'^'^''U{t) | Ve^en)^ =F {ipevenMt) l^Jeven) p ■ (^19) 

with 



exp ^— i j dtHoF{t)^ 



(B20) 



U{t) = 
We denote 

iloFit) = e'^^^y^+^-^*noFe-'^'^'y'+'^-^' = ^ (o{t) \n + l)p p {n\ e''^* + d\t) \n)p p{n+l\ e"''^*) . (B21) 

n 

with 6{t) = e^'^'y'*de-''^'y'\ 

Now we are ready to expand the interaction Hamiltonian in powers of T-LoFit)- Up to second order we encounter 
the terms: 

U{t) - 1 = - /o* dt2 Jo' dti (d{t2) \n + l)p e*("+iV*= + d^{t2) \n - l)p e*("-i)"*2^ {p (n| n)Fe-'""(*^-*i)) 

(dHh)F {n + 1| e-^("+i)'^*i + 6{ti)p (n - 1| e-'("+i)'^*i) . ^^^^^ 

Note that we split the compound: \n)p p (m| e-*'^(™-")* = \n)p e*'^"* -p {m\ e"*'^'"*. The operators 6, could at 
this point be construed as first quantized operators, which change the state of a particle interacting with the radiation 
field. 

Further progress is made by projecting on the initial, intermediate, and final states of the system described by 
Tigys- Let us denote P„ to be a projector of the system's state on the subspace of energy e^r- For now, we maintain the 
generality of the discussion, although eventually, wc will restrict ourselves to Hgys = H{\i.) which is a 4x 4 Hamiltonian 
describing two 2d subspaces with energies ±e(k); at that point it will be simple to use a = ±1 to indicate the valence 
vs. conduction subspaces. Armed with this notation we can write: 

U(t)-1 = - J2 Jo dt2 J^' dhP^^ (6(t2)|n+l)^e*("+i)^*^+Ot(i2)|n-l)^e*("-i)"*^)p<,,e-'"'^(*^-*i) 

(d\ti)F {n + 1| e-«^("+i)*i + d{ti)F {n - 1| e-^'^("-i)*i) 

(B23) 

This allows us to resolve the time dependence of the operators: 

L{{t) - 1 = - E Jo dt2 /o' dtiP^.e-'^'"^-^"^^*^ (o \n + l)p e^("+i)'^*= +d^\n- l)p e*("-i)'^*2] P^^ 

<Tl,a2,<T3 ^ ^ (B24) 

id'^p {n + 1| e-'("+i)"*i +dF{n- 1| e-*^"-!)'^*! j p^^e-*('=''i-'=''2)*ie-"'^(*2-*i) 

The expressions compactify by defining two indices /xi^2 = ±1, and denoting 0^+^^ = and 0^~^^ = O, and dropping 
the subscript F, 

AU{t)-i = - E E fdt2fj^dhP,,d(-^^) 

<7i,<72,0-3 jUi,/:i2=±l (B25) 
|j^_|.^2) P„ O^**^^ (n + /ill Per e~'^*''l~*''2^*^~'^""'"''^'"*^e~'^^''2~^<'3)*2+»("+A'2)wt2g-mw(t2-tl) 

By moving to average time, t = and time difference, t- = t2 — ii , as well as integrating over t_ (while assuming 

that t is large and ignoring boundary terms for t-, we get: 

U{t)-l = - E dtP,,0^->^^^\n + H2)P..0^^^Hn + ^i\Pa. ^^^^ (B26) 

o-i, 0-2, 0-3 jui,/:i2=±l 2 "1"'^ 2 ^<^2 

The time dependence on t is simply the interaction-representation time dependence. Therefore, by going back to the 
Schrdinger representation, we are able to get rid of the remaining time dependence in the expression, and we readily 
extract the effective second-order contributions to the effective I-Lf'- 

n'-U= E E El- + /^^)..<- + Mr|gg— (B27) 

o-i, 0-2, era /ni,/t2=±l n 2 2 ^<^2 
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whereby now 



ftp — tts 



(B28) 



Again we note that \il^even) p is an eigenstate with eigenvalue 1 of ^ |n + H2)f f ('^ + ^oi 1x1^2 — ±1, which allows 

n 

the mapping back to the original system. As before, the way to go back to the language of the original time dependent 
problem is to evaluate Eq. (B18), which yields 



U2{t) 



exp 


h/ 




. Jo 



dt 



Jjeff 



it)) 



(B29) 



with 



h: 



eff 
2-ph 



it) 



<7l, <T2, <T3 fJ.i^lJ,2=:tl n 



(B30) 



The form of Eq. (B30) is clearly in accord with degenerate perturbation theory. The reason for the putative 
degeneracy is the fact that we consider the energy of the F states representing the photons together with the energy 
of the system. A resonance, therefore, translates to a degeneracy in this language. It is interesting to note that Eq. 
(B27) generalizes degenerate perturbation theory to the case of near degeneracy. The energy denominator is actually 
the difference between the average of the initial and final energies, and the intermediate energy. 



2. Application to the 3d FTI construction 

Applying the formalism above to the 3d FTI construction is quite straightforward. We let O = = = ^V-'y, 
and Pa = 5 ^1 + "'^^p) with -ff(k) = D ■ ^. We separate to two cases: (1) ai = 0-3, fii = fj,2, and (2)cri = —as, 
Ml = -/"2 = -0-1. 

a. Case 1 - diagonal elements. The case of diagonal elements can be treated for both the valence and conduction 
band simultaneously, since terms that do not excite between the bands are time independent, and therefore the f terms 
factor out. Therefore: 

CTl,CT2=±l, ([i=±l ^ ' 

Clearly a2 = —cri, otherwise the denominator makes the sum vanish. Thus: 
Elementary algebra of the Dirac matrices reduces this expression to: 



(v- ^d-vd)\ 

K{U = f {n\ 4,2_^2 ^ 2^ • 7 (B33) 

where we also recognize Vj_f) = V — -^D ■ VD. Note that a term corresponding to Eq. (B33) also arises in the 
treatment involving the two rotating wave transformations. Consider the Hamiltonian H2{t), Eq. (AlO), evaluated at 
values of k for which 2w = e(k), in the rotating wave approximation. For these k values, the terms s'"^\'k)P^\]s.) + 

(e!^^(k) + w)Pl^^(k) do not vanish, but give, to order \Vi\^ a term corresponding to Eq.(B33). 

b. Case 2 - interhand elements. The interband elements are to some extent more complicated, since we need to 
consider excitation and relaxation separately. We consider a term connecting the initial state cti with 0-3 = For 
this process to be viable, we must have = —^2 = — fi- Therefore, the specific term is: 

K{U-.,=\n + i,2)pAn + liA E (B34) 
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Using the form of Pct2 we can carry out the sum, and obtain: 



Once again, elementary manipulations of the Dirac matrices yields: 



K{U-., = |„ + ^2)^ ^ (n + Ml I (B35) 



K%-a. =\\n + M2)^ ^ (n + Mil ( -^^±6 • 7 - p . 7, ^ • 7j 1 . (B36) 



This term seems indeed complicated. A simphfication occurs, however, when we consider the RWA directly applied 
with a 2oj energy. The effective Hamiltonian connecting the two bands (that arises from V) is (the RWA eliminates 
the time dependence, and hence the F's): 



"H^iUan. - KH-, + n'^J'Ll (B37) 



This term joins the direct two-photon process that arises from the term appearing in due to the mass curvature. 
This term also coincides to second order in V, with the corresponding effect on the radiation-induced gap within the 
two consecutive rotating wave transformations. Note that at resonance Ck = 2Aek, and thus Eq. (B37) agrees with 
Eq. (A19). 
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